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(57) ABSTRACT 

An airdata estimation and evaluation system and method, 
including a stable algorithm for estimating airdata from 
nonintrusive surface pressure measurements. The airdata 
estimation and evaluation system is preferably implemented 
in a flush airdata sensing (FADS) system. The system and 
method of the present invention take a flow model equation 
and transform it into a triples formulation equation. The 
triples formulation equation eliminates the pressure related 
states from the flow model equation by strategically taking 
the differences of three surface pressures, known as triples. 
This triples formulation equation is then used to accurately 
estimate and compute vital airdata from nonintrusive surface 
pressure measurements. 

17 Claims, 8 Drawing Sheets 
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STABLE ALGORITHM FOR ESTIMATING 
AIRDATA FROM FLUSH SURFACE 
PRESSURE MEASUREMENTS 

ORIGIN OF THE INVENTION 5 

The invention described herein was made in the perfor- 
mance of work under a NASA contract, and is subject to the 
provisions of Public Law 96-5 17 (35 U.S.C. § 202) in which 
the Contractor has elected not to retain title. 
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BACKGROUND OF THE INVENTION 

1. Field of the Invention The present invention relates in 
general to flow instrumentation and more particularly to an 
airdata estimation and evaluation system and method for 
estimating and evaluating airdata from noninstrusive surface 
pressure measurements. 

2. Background Art 

Airdata are critical information for any type of flight 
vehicle. Airdata are those parameters, characteristics, prop- -0 
erties and quantities derived from the air surrounding the 
flight vehicle. Accurate airdata are absolutely necessary for 
many purposes and applications and help to ensure efficient 
and safe flight. 

For example, a pilot cannot safely fly an aircraft without 25 
knowing its airspeed and pressure altitude. In civil aviation, 
air traffic controllers determine the separation distances 
between aircraft based on accurate knowledge of the air- 
craft's pressure altitude. If the airdata is inaccurate, the 
aircraft’s true pressure altitude may differ greatly from its 30 
reported pressure altitude, and disaster may result. 

Numerous other aircraft systems also need accurate air- 
data. Autopilot, engine controls, cockpit and cabin environ- 
mental controls, weapons delivery, navigation and air traffic 
control are just a few examples of aircraft systems that 35 
depend on accurate airdata. 

Airdata are most accurate when the undisturbed air away 
from the flight vehicle is measured. This is because the 
presence of the flight vehicle in the air can cause measure- 
ment errors. Moreover, as the velocity of the flight vehicle 
increases, compressibility and shock waves further disrupt 
the surrounding air. 

Traditionally, in order to mitigate these errors, airdata 
measurements have been performed using booms that 45 
extend outward from the flight vehicle and measure the 
undisturbed air. Typically, these intrusive booms contain a 
pressure-measuring instrument that operates by stagnating 
the freestream flow and measuring the difference between 
static and impact pressure. These booms are usually located 50 
at the nose of the flight vehicle. Although these booms are 
excellent at making steady measurements at low-to- 
moderate angles of attack, the booms are sensitive to 
vibration, alignment error and are easily damaged in flight or 
on the ground. 55 

There are also several applications where an intrusive 
boom is undesirable. For example, on hypersonic aircraft an 
intrusive boom may either be tom off by the forces generated 
in flight or melted by the high temperatures involved in 
high-speed flight. On stealth aircraft, where a minimum 60 
radar cross section is required, an intrusive boom is unde- 
sirable because it makes the aircraft highly visible to radar. 
Therefore, intrusive booms are generally unsuitable for 
airdata measurements in more advanced types of flight 
vehicles. 65 

One technique developed to overcome these difficulties 
with the conventional intrusive boom was a flush airdata 


2 

sensing (FADS) system. Unlike the intrusive boom that 
directly measures airdata from the freestream air, the FADS 
system infers airdata from nonintrusive surface pressure 
measurements. 

The fundamental concept of the FADS system is that it 
does not require probing of the freestream flow in order to 
compute airdata. Instead, airdata can be estimated and 
evaluated from flush surface pressure measurements. These 
pressure measurements are then related to airdata using 
aerodynamic model equations that mathematically describe 
the flow. Solving the aerodynamic model equations math- 
ematically extracts the airdata given the surface pressure 
measurements. 

One early use of a FADS system w as in the 1960s by the 
X-15 high -altitude, hypersonic aircraft. The X-15 FADS 
system used a hemispherical nose that was actively steered 
into the local relative wind to make airdata measurements. 
One drawback to this system, however, was that the actively 
steered nose was an extremely complicated mechanical 
design. Consequently, the steered-nose FADS system was 
abandoned with the termination of the X-15 program. 

A modem FADS system was developed as an experiment 
on the space shuttle. The shuttle entry airdata system 
(SEADS) was a k *ride -along” experiment aboard the shuttle 
but never was used to actually provide airdata to the shuttle 
pilot. The SEADS used a matrix of fixed static-pressure 
measurements and required no mechanical steering of the 
nose. 

The problem, however, with both of these systems is that 
both are quite inaccurate. This is because the two systems 
merely focused on the feasibility of the FADS concept and 
made no attempt to derive algorithms for more accurate 
estimation and evaluation of airdata from the surface pres- 
sure measurements. 

The next advance in FADS technology was a real-time, 
high angle-of-attack FADS (HI-FADS) system. The 
HI-FADS design, as with the earlier X-15 FADS system and 
the SEADS, used a matrix of flush static-pressure orifices 
arranged on the nose of the flight vehicle. Moreover, the 
HI-FADS system did attempt to derive algorithms for esti- 
mating and evaluating airdata from surface pressure mea- 
surements. 

One problem, however, with the HI-FADS design is the 
complexity of the algorithm and the associated increased 
costs. In particular, the system incorporates all surface 
pressure measurements simultaneously into an overdeter- 
mined estimation algorithm. Furthermore, airdata param- 
eters are inferred using nonlinear regression methods. 
However, the use of these nonlinear regression cause the 
estimation algorithm to be highly unstable and complex. In 
addition, the cost of maintaining the software implementing 
this algorithm is quite high. 

In order to solve these aforementioned problems with 
these prior art FADS systems, the present invention was 
developed to provide a reliable and accurate FADS system 
for use in the X-33 flight vehicle program. The X-33 is 
53-percent scale model of the Single Stage to Orbit (SSTO) 
Reusable Launch Vehicle (RLV) rocket system. The primary 
goal of the SSTO program is to radically reduce the cost of 
access to space, and the X-33 is a demonstrator of this 
advanced technology. The X-33 is designed to achieve a 
peak altitude near 300,000 feet and speeds of greater than 
Mach 12. After atmospheric re-entry, the X-33 returns to 
Earth in an unpowered horizontal landing. 

Because the X-33 will perform an unpowered landing (i.e. 
there is only a single landing attempt), highly accurate 
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determination of the airdata such as dynamic pressure, 
angle-of-attack, and surface winds is critical to insure that 
the target runway can be reached under a wide variety of 
atmospheric conditions. Furthermore, direct feedback of 
angle-of-attack and angle-of-sideslip may reduce the gust- 
load on the vehicle airframe during the ascent phase of the 
flight. 

The determination was made early in the X-33 program 
that the full airdata state including Mach number, angle-of- 
attack, angle-of-sideslip, dynamic pressure, airspeed and 
altitude would be a critical requirement for both the RLV and 
the X-33 and that existing FADS algorithms could not 
adequately fulfill this requirement. 

Specifically, there are several problems with existing 
FADS algorithms. One problem is that the nonlinear regres- 
sion method used in the estimation algorithm tends to be 
highly unstable. In fact, the software implementing this 
algorithm requires numerous special software patches in an 
attempt to maintain stability. 

This highly unstable algorithm leads to the problem of 
complexity. Specifically, the special software patches 
required to make the algorithm stable also require a complex 
maze of software designed to maintain stability. Moreover, 
the cost of establishing and maintaining this software is 
extremely high. 

Another problem with existing FADS systems is that the 
unstable estimation algorithm can lead to numerous reliabil- 
ity problems. This, in turn, has detrimental consequences for 
the flight vehicle that relies on accurate airdata. 

Therefore, what is needed is a FADS system that has a 
highly stable estimation and evaluation algorithm that can 
easily be implemented with a minimum of cost and main- 
tenance. Moreover, what is further needed is a FADS system 
having an estimation and evaluation algorithm that is 
extremely reliable and provides highly accurate airdata for 
use by the various systems of the flight vehicle. 

Whatever the merits of the above-mentioned FADS sys- 
tems and methods, they do not achieve the benefits of the 
present invention. 


NOMENCLATURE 


A 

Angle -of- Attack Triples Algorithm Coefficient 

a 

FADS geometry Coefficient 

A' 

Angle -of- SidesUp Triples Algorithm Coefficient 

Ao(M^) 

Zero’ th Order Angle of Attack Calibration Coefficient 

A,(M^) 

First Order Angie of Attack Calibration Coefficient 

A4M_) 

Second Order Angle of Attack Calibration Coefficient 

A,(MJ 

Third Order Angle of Attack Calibration Coefficient 

B 

Angle-of- Sideslip Triples Algorithm Coefficient 

B 0 <M_) 

Zero' th Order Angle of Sideslip Calibration Coefficient 

B,(M_) 

First Order Angle of Sideslip Calibration Coefficient 

B4M_) 

Second Order Angle of Sideslip Calibration Coefficient 

B,<M_> 

Third Order Angle of Sideslip Calibration Coefficient 

b 

FADS Geometry Dummy Variable 

B’ 

Angle-of- SidesUp Triples Algorithm Dummy Variable 

C 

Angle-of- Sideslip Triples Algorithm Dummy Variable 


Subsonic Mach Number Function 

fi 

Supersonic Mach Number Function 

M 

Estimation Algorithm Geometry Matrix 

M- 

Free stream Mach number 

P- 

Freesueam Static Pressure, psf. 

P. 

Pressure Value, i'th FADS port, psf. 

Q 

Pressure Weighting Matrix 

4 

Weights on Measured Pressures 

qc 

Impact pressure, psf. 

X 

State Vector 

X : 

Mean Squared Fit Error, psf* 
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Angle of attack, generic, degrees 

Local Angle of attack sensed by FADS, degrees 

Frerstream Angle of attack, degrees 

Angle of sidesUp. generic, degrees 

Local Angle of sideslip sensed by FADS, degrees 

Frerstream Angle of sideslip, degrees 

Quadratic Solution 1 for Angle of sideslip, degrees 

Quadratic Solution 2 for Angle of sideslip, degrees 

Angle of attack correction term, degrees 

Angle of sideslip correction term, degrees 

Pressure Difference for Triples Algorithm, psf 

Ratio of Specific Reals 

Position Error Calibration Parameter 

Variation of € with Mach number at a e , * 0 

First Order Coefficient, Fit of t with Angle of Attack 

Second Order Coefficient, Fit of « with Angle of Attack 

First Order Coefficient, Fit of « with Angle of Sideslip 

Second Order Coefficient. Fit of « with Angle of Sideslip 

Cone angle of FADS port, degrees 

Clock angle of FADS port, degrees 

Local flow incidence angle, degrees 

Standard deviation 

Dummy Variable Used in Angle of Sideslip Root 
Analysis 

Chi Squared Distribution 

FADS Geometry Term, cos'9 + « sin^ 
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Superscripts, Subscripts, and Mathematical Operators 


30 


35 


40 


45 


i 

Port index 

j 

Port index 

0> 

Iteration Index 

k 

Port index 


FADS Estimate 

Det(.J 

Matrix Determinant 

ir‘ 

Matrix Inverse 

i.r 

Matrix Transpose 

VM, 

Gradient with Respect to Vector x 

dl) 

Partial Derivative with Respect to x 

Ox 


Acronyms 


A/D 

Analog to Digital Conversion 

FADS 

Flush Airdata Sensing 

RLV 

Reusable Launch Vehicle 


SUMMARY OF THE INVENTION 

50 

To overcome the limitations in the prior art as described 
above and other limitations that will become apparent upon 
reading and understanding the present specification, the 
present invention is embodied in an airdata estimation and 
55 evaluation system and method for estimating and evaluating 
airdata from nonintrusive surface pressure measurements. 

The airdata estimation and evaluation system of the 
present invention includes a triples formulation module for 
eliminating the pressure related states from the flow model 
60 equation, an angle of attack module for computing the angle 
of attack, an angle of sideslip module for computing the 
angle of sideslip and an airdata module for estimating and 
evaluating other airdata. Furthermore, the present invention 
includes a method for estimating and evaluating airdata 
65 using this apparatus. The system and method of the present 
invention determine airdata in an economical, reliable and 
highly stable manner. 
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In a preferred embodiment, the airdata estimation and 
evaluation system is utilized in a FADS system on a vehicle 
within a flowfield. Preferably, the FADS system includes 
nonintrusive flush-mounted surface pressure sensors for 
providing pressure measurements and the airdata estimation 
and evaluation system for computing the airdata parameters. 
Moreover, the computed airdata is received by an airdata 
processor that determines whether, based upon the airdata 
received, any changes or adjustments to the vehicle are 
necessary. These adjustments may include changes in veloc- 
ity and attitude. If the airdata processor determines that 
changes are necessary, these changes are effected by a 
control system in response to the airdata processor’s com- 
mands. The control system may make effect changes by 
physically moving the vehicle’s control surfaces or adjusting 
the power delivered by the powerplant. 

The airdata estimation and evaluation system includes a 
triples estimation module, for evaluating a triples formula- 
tion equation using the surface pressure measurements and 
an angle of attack (AOA) module. The AOA module is used 
for computing an angle of attack using the triples formula- 
tion equation. Furthermore, the triples estimation module 
includes an angle of sideslip (AOS) module, for computing 
the angle of sideslip using the computed angle of attack and 
the triples formulation equation. Further, the system 
includes an airdata module for estimating and evaluating 
other airdata using the computed angle of attack and the 
computed angle of sideslip. 

The airdata estimation and evaluation system of the 
present invention may also be implemented as a method as 
a dedicated logic circuit or a field programmable gate array. 
Preferably, the airdata estimation and evaluation system 
contains a microprocessor and memory for storing the 
instructions coded with the method of the present invention. 

The method of the present invention includes a triples 
formulation function, an angle of attack function, an angle of 
sideslip function and an airdata function. The triples formu- 
lation function transforms the given flow model equation 
into a triples formulation equation. 

The angle of attack function computes an angle of attack 
from the triples formulation equation. In a similar manner, 
the angle of sideslip function computes the angle of sideslip. 
The airdata function estimates and evaluates airdata using 
the flow model equation, the computed angle of attack and 
the computed angle of sideslip. 

The present invention may also be implemented in a 
computer-readable medium whose contents cause a com- 
puter to function as an airdata estimation and evaluation 
system as described above. Moreover, the present invention 
may also include a computer-readable memory device 
encoded with a data structure for estimating and evaluating 
airdata as described above. This airdata used by the airdata 
processor to determine the necessary corrections to the 
vehicle in the flowfield. 

Other aspects and advantages of the present invention as 
well as a more complete understanding thereof will become 
apparent from the following detailed description, taken in 
conjunction with the accompanying draw ings, illustrating by 
way of example the principles of the invention. Moreover, it 
is intended that the scope of the invention be limited by the 
claims and not the preceding summary or the following 
detailed description. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Referring now to the drawings in which like reference 
numbers represent corresponding parts throughout: 


6 

FIG. 1 is an overview of the flush airdata sensor (FADS) 
system of the present invention as implemented on an X-33 
flight vehicle. 

FIG. 2 is a front view of the X-33 flight vehicle of FIG. 
5 1 illustrating the preferred configuration of the surface 
pressure sensor matrix on the nose of the vehicle. 

FIG. 3 is a side view of FIG. 2 further showing the 
configuration of the six sensors of the surface pressure 
sensor matrix. 

FIG. 4 a perspective view of the nose of the X-33 flight 
vehicle in FIGS. 2 and 3 illustrating the coordinate angles 
for the surface pressure sensor matrix. 

FIG. 5 is an overview structural block diagram of the 
15 airdata estimation and evaluation system shown in FIG. 1. 

FIG. 6 is a detailed functional flow diagram of the triples 
formulation function of the airdata estimation and evaluation 
system of FIGS. 1 and 5. 

FIG. 7 is a detailed functional flow diagram of the angle 
20 of attack function of the airdata estimation and evaluation 
system of FIGS. 1 and 5. 

FIG. 8 is a detailed functional flow diagram of the angle 
of sideslip function of the airdata estimation and evaluation 
system of FIGS. 1 and 5. 

~ 5 FIG. 9 is a detailed functional flow diagram of the airdata 
function of the airdata estimation and evaluation system of 
FIGS. 1 and 5. 

DETAILED DESCRIPTION OF THE 
30 PREFERRED EMBODIMENT 

In the following detailed description of the preferred 
embodiment, reference is made to the accompanying draw- 
ings in which is shown by way of illustration a specific 
35 embodiment whereby the invention may be practiced. It is to 
be understood that other embodiments may be utilized and 
structural changes may be made without departing from the 
scope of the present invention. 

I. Introduction 

40 As shown in the drawings for purposes of illustration, the 
present invention is embodied in an apparatus that includes 
an airdata estimation and evaluation system. The invention 
may also be implemented as a method for estimating and 
evaluating airdata using this system. The system and method 
45 of the present invention provides a highly accurate, 
extremely reliable and relatively simple way for airdata to be 
extract from nonintrusive surface pressure measurements. 
Existing FADS systems include highly unstable algorithms, 
extremely complicated software implementing these algo- 
50 rithms and are inaccurate and unreliable. 

The present invention includes an airdata estimation and 
evaluation system for estimating and evaluating airdata from 
surface pressure measurements. This system includes a 
triples formulation module for transforming the flow model 
55 equations into a triples formulation equation in order to 
eliminate certain variable from the flow model equation. The 
system also includes an angle of attack module for comput- 
ing an angle of attack from the triples formulation equation 
and an angle of sideslip module for computing the angle of 
60 sideslip from the triples formulation equation. Furthermore, 
the system includes an airdata module for computing the 
remaining airdata quantities using the angle of attack and the 
angle of sideslip. 

The method of the present invention may be implemented 
65 in the above system and includes a way to estimate and 
evaluate airdata from nonintrusive surface pressure mea- 
surements using the aforementioned apparatus. 
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It should be noted that, although throughout specification 
reference is made the air as the medium, nearly any fluid or 
gas capable of having pressure measured could be used in 
place of air. In this situation, the airdata would become 
airdata, since air would no longer be the medium. For 
example, the present invention may be utilized in a subma- 
rine moving through water. 

The airdata estimation and evaluation system of the 
present invention is capable of accurately computing airdata. 
In particular, the surface pressure data is strategically chosen 
such that certain pressure dependent terms drop out of the 
flow model equation. This permits the flow model equation 
to be solved without the use of nonlinear regression methods 
and results in a more accurate and a more stable algorithm. 
Moreover, the code required to implement and maintain the 
system of the present invention is more economical and less 
complex than prior art systems. 

II. Structural Overview of the Invention 

FIG. 1 is an overview of the flush airdata sensing (FADS) 
system 100 suitable for utilizing the new airdata estimation 
and evaluation system and method of the present invention. 
The FADS system 100 of the present invention is shown 
implemented on a flight vehicle 110. In the illustration, the 
flight vehicle 110 is the X-33 test vehicle as described above. 

The FADS system 100 includes a surface pressure sensor 
matrix 120 located at the nose area of the flight vehicle 110. 
Preferably, the number of sensors in the surface pressure 
sensor matrix 120 is six. This number was selected as a 
compromise between the need to accurately measure the 
flow conditions at the nose and the cost of locating the 
sensors on the flight vehicle 110. Because, in the case of a 
flight vehicle, there are four airdata states and a calibration 
parameter to be estimated, at least five independent pressure 
measurements must be available to derive the entire airdata 
state. Using five sensors in the surface pressure sensor 
matrix 120 to estimate the airdata is equivalent to a higher 
order spline fit and results in an estimation and evaluation 
algorithm that is sensitive to noise in the measured pres- 
sures. However, if an additional sixth sensor is provided, the 
result is to mitigate the noise sensitivity, increase redun- 
dancy options, and results in a FADS system 100 that gives 
overall superior performance. 

FIG. 2 is a front view of the X-33 flight vehicle 110 of 
FIG. 1 illustrating the preferred configuration of the surface 
pressure sensor matrix 120 on the nose 205 of the vehicle 
110. Moreover, FIG. 3 is a side view of the nose 205 in FIG. 
2 further showing the location of the six sensors of the 
surface pressure sensor matrix 120. 

Referring to FIGS. 2 and 3, five sensors are located on the 
nose 205 and a sixth sensor 210 is located slightly under the 
nose 205. Four of the sensors, including the sixth sensor 210, 
are located along the meridian 220. As further explained 
below, this layout permits the calculations for angle of attack 
to be decoupled from the calculations for angle of sideslip 
using the method and apparatus of the present invention. In 
addition, this sensor layout is designed to give good sensi- 
tivity for local angles of attack varying from -20° to 45°, 
and angles of sideslip up to ±20°. 

FIG. 4 is a perspective view of the nose 205 of the X-33 
flight vehicle 110 shown in FIGS. 2 and 3 illustrating the 
coordinate angles for the surface pressure sensor matrix 120. 
The cone angle 410 (X) is the total angle the normal to the 
surface 420 makes with respect to the longitudinal axis 430 
(x axis) of the flight vehicle 110. The clock angle 440 (<J» is 
the clockwise angle looking aft around the axis of symmetry 
450 starting at the bottom of the flight vehicle 110. Using 
these coordinate angles, the nominal clock and cone angles 


8 

of the six sensors in the surface pressure sensor matrix 120 
are tabulated in Table 1 : 


TABLE I 


FADS Purls CI<Kk ami Cone Annies 
Port # 0, 

X, 

1 

181) 

20 

-> 

270 

20 

3 

0 

0 

4 

90 

20 

5 

0 

20 

6 

0 

45 


15 Referring again to FIG. 1, the above -described surface 
pressure sensor matrix 120 is connected to a pressure 
distribution interlace 130. Typically, this interface 130 is a 
device that transforms analog measurements into digital 
data, such as an analog to digital (A/D) converter. However, 
20 this example is only one of several possible configurations 
that one skilled in the art would recognize. 

Through the interface 130, the surface pressure sensor 
matrix 120 is connected to the airdata estimation and evalu- 
ation system 140. This system 140 inputs measurements 
25 from the surface pressure sensor matrix 120 and estimates 
and evaluates airdata. 

This computed airdata is sent to an airdata processor 145 
that determines whether a physical change in the parameters 
the flight vehicle 110 are necessary. If so, then the control 
30 system 150 instructs the attitude controller 160, the velocity 
controller 170 or other vehicle system to effect the change. 
For example, the airdata processor may instruct the control 
system 150 to change the flight vehicle's 110 orientation, 
velocity or other physical parameter requires adjustment. 

35 In response to the decision of the airdata processor 145, 
the control system 150 sends control signals to one or more 
of the attitude controller 160 and the velocity controller 170. 
Both of these controllers may be used to physically change 
the control surfaces, engine thrust or other aspect of the 
40 flight vehicle's 110 flight parameters. Moreover, the airdata 
processor 145 may send the computed airdata to the instru- 
mentation system and navigational system of the flight 
vehicle 110. Several other uses of the computed airdata are 
possible, as can be appreciated by those of ordinary skill in 
45 the art. 

1 . Airdata Estimation and Evaluation System 

FIG. 5 is an overview structural block diagram of the 
airdata estimation and evaluation system 140 shown in FIG. 
1. The system 140 includes a triples formulation module 
50 510, an angle of attack (AO A) module 520, an angle of 
sideslip (AOS) module 530 and an airdata module 540. As 
described above, the inputs to the system 140 are the 
pressure distribution interface 130 and the output of the 
system 140 is the airdata processor 145. 

55 The purpose of the airdata estimation and evaluation 
system 140 is to evaluate and estimated airdata given the 
surface pressure sensor matrix 120. This estimation and 
evaluation is accomplished by extracting airdata from the 
nonintrusive surface pressure measurements using an aero- 
60 dynamic or flow model equation. In order to be useful, this 
flow model equation must capture the salient features of the 
flow and be valid over a large Mach number range. 
Moreover, the flow model equation must be simple enough 
to be inverted in real-time so that the airdata can be 
65 extracted. 

To solve the problem of describing a complex flow 
scenario with a simple model, the flow model equation of the 
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present invention was derived as a splice of the closed-form 
potential How solution for a blunt body, applicable at low 
speeds. Furthermore, a modified Newtonian flow model, 
applicable at hypersonic speeds, was spliced with the poten- 
tial flow solution. Both the potential flow and modified 5 
Newtonian flow describe the measured pressure coefficient 
in terms of the local surface incidence angle. 

In order to blend these two solutions over a targe range of 
Mach numbers, a calibration parameter E was used. This 
parameter must be empirically calibrated to allow for the 10 
effects of flow compression, body shape, and other system- 
atic effects such as shock wave compression or Prandtl- 
Meyer expansion on the forebody. However, once s is 
calculated it can be use for that particular flight vehicle. 

The flow model equation for the present invention is: 15 

sin : (8pF-^~ 0) 

Since equation (1) is derived from the potential and 
Newtonian flow around a blunt body, it is most valid near the 
vehicle stagnation point. Accordingly, the preferred location 
for the surface pressure sensor matrix 120 is at the nose of 
the flight vehicle 110. 

In equation { 1 ), 0, is the flow incidence angle between the 
surface normal at the \ th port and the velocity vector. The 
flow incidence angle is related to the local (or effective) 
angle-of-attack a e , and angle-of-sideslip (3,. by the equation: 

cos(0 / >=cos<a ( .)co&(^)cos<k ; >-Kin(^)sin(4> I )sm(X / H 

sinio^Kosi^Jcos^sinik,) (2) 

30 

Referring to FIG. 4 and equation (2), the cone angle X is 
the total angle the normal surface makes with respect to the 
longitudinal axis of the nosecap. The clock angle <|> is the 
clockwise angle looking aft around the axis of symmetry 
starting at the bottom of the fuselage. The remaining param- 35 
eters in equation (1) are impact pressure qc and the 
freestream static pressure P^. Using these four basic 
parameters, a r , (3,, q c , and P^, most other airdata quantities 
of interest may be directly calculated. 

In addition to the calibration for e, the local flow inci- 40 
dence angles, a,, and (3,, must be related by calibrations to 
the true freestream flow angles. These flow-angle calibra- 
tions account for such additional systematic effects as bow 
shock flow deflection and body-induced upwash and side- 
wash. A preliminary set of calibrations for the X-33 flight 45 
vehicle has been derived using wind tunnel data. In a similar 
manner, these calibrations may be empirically determined 
for other types of flight vehicles, 
a. Triples Formulation Module 

The flow model equation of equations (1) and (2) is 50 
inherently non-linear in terms of the state parameters and 
can not be directly inverted. As discussed in the Background 
section, a non-linear regression method in which the equa- 
tions are recursively linearized and inverted using iterative 
least-squares has been used with marginal success. A major 55 
drawback, however, is that the algorithm is highly unstable 
and thus requires complex software patches just to maintain 
stability. 

The present invention includes a new system and method 
to solve equations 1 and 2 to thereby extract airdata. One 60 
way this is accomplished is by using the airdata estimation 
and evaluation system 140. At the heart of the airdata 
estimation and evaluation system 140 is the triples formu- 
lation module 510. 

Referring again to FIG. 5, the triples formulation module 65 
510 is a part of the airdata estimation and evaluation system 
140. The triples formulation module 510 takes the flow 


10 

model equation (equation ( 1 )) and transforms it into a triples 
formulation equation. This is done by taking a strategic 
combination of three surface pressure readings from the 
surface pressure sensor matrix 120 such that the flow 
angularity states are decoupled from the static and impact 
pressure states and the aerodynamic calibration. 

Significantly, taking the differences of this strategic selec- 
tion of three surface pressures (referred to as “triples”) 
eliminates the pressure related states, q c , P w , and e from 
equation ( I) (the flow model equation). This is a significant 
and novel discovery permitting a more stable and direct 
solution of equation (1) than achieved by existing algo- 
rithms. The mathematics of the triples formulation module 
510 as well as of all the other modules contained in the 
airdata estimation and evaluation system 140 is detailed 
below in the “Mathematical Description of the Invention” 
section. 

After the flow model equation has been transformed, the 
resulting triples formulation equation is: 

r.* cos*0y+r y7 cos^+r*, cos^o^o (3) 

where, 

fariPrPk)* Tj&PrPi)- V k^Pk~Pj> < 4 > 

and p ( , p jy and p k are the pressures used in the data triple, 
b. Angle of Attack (AOA) Module 

The triples formulation equation is sent to the AOA 
module 520, for computation of the angle of attack a. This 
is calculated using the triples formulation equation and the 
pressure quantities supplied by the surface pressure sensor 
matrix 120 through the pressure distribution interface 130. 

In equation (3), the local angle of attack can be decoupled 
from the local angle of sideslip by using only pressure ports 
aligned along the vertical meridian 220 where <J>=0 o or 
<(>= 1 80°. In this geometry arrangement terms related to angle 
of sideslip fall out of equation (3). Thus, for |a|^45° the 
solution for a is: 



where 

A=r* siirXy+T), siirX^+r*, sin 2 X 

cos sin Xj cos cos sin A* cos ~^ k +T k> cos 4, sin 
A., cos A., (6) 

When |a|>45 0 the correct solution is given by the comple- 
ment of equation 5: 



For the surface pressure sensor distribution 120 there are 
four possible independent combinations of a-triples on the 
vertical meridian. The clock and cone angles corresponding 
to these triples combinations are presented in Table II: 
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TABLE II 


Clock ami Cone Angles 

of the An*jle-of- Attack Triples 

K- 

Triple # 

<fc 

Ay. 

4>, 

K 

♦k- 

, 

0 

0 

0 

20 

180 

20 


0 

0 

0 

20 

0 

45 


0 

0 

180 

20 

0 

45 

4 

0 

20 

180 

230 

0 

25 


In the preferred embodiment, -the output angle of attack is 
the mean of the angle of attack values computed using the 
four individual triples. This averaging procedure provides a 
measure of noise rejection for the airdata estimation and 
evaluation system 140. Clearly, if one of the sensors in the 
vertical meridian 220 of the surface pressure sensor matrix 
120 is deemed unusable and is weighted out of the 
algorithm, then only one valid triple remains for computing 
the output angle of attack, 
c. Angle of Sideslip (AOS) Module 

After the angle of attack a has been estimated, the angle 
of sideslip p is evaluated by the AOS module 530. This 
module uses the triples formulation equation from the triples 
formulation module 510, the output angle of attack from the 
AOA module 520 and pressure quantities from the surface 
pressure sensor matrix 120 to calculate the angle of sideslip. 

The AOS module 530 computes the angle of sideslip by 
using a combination of the available pressure sensors in the 
surface pressure sensor matrix 120, other than the obvious 
set in which all three sensors lie on the vertical meridian. 
The result is a quadratic equation in tan p; 

A' tan 2 $+2ff tan f^C'=0 (8) 

where 

B = ( r lk u l b^T jl a k b l ^Y i t /jAl 

C'{r^/+r^, 2 +r^r} (9) 
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Consequently, the solutions for p_ and q ( must be extracted 
iteratively. These iterative computations are performed by 
the airdata module 540. The airdata module 540 accepts 
input from both the AOA module 520 and the AOS module 
s 540. 

Defining the Matrices: 


to 


COs'ltf] ) +■ \ 

) J 

cos : <0 ft ) + \ 

^ ;, sirr(0 6 > J 




(in 


the flow model equation (equation (1)) can be used to 
develop an iterative estimator of the form: 


<ic 


P\ 

Poo 

t j* l ) 

P6 


The subscript (j) refers the result of the j ,f/l iteration. The 
notation M (y) refers to the matrix of equation 12 with € being 
evaluated using the values for Mach number resulting from 
the previous iteration. The q ; terms are weights which have 
a nominal value of 1 .0. Setting the value of q, to zero weights 
the i th pressure reading out of the algorithm. 

Given q c and p„, the Mach number can be computed 
using normal one-dimensional fluid mechanics relation- 
ships. For subsonic speeds, the Mach number can be com- 
puted using isentropic flow- laws, where for y= 1 .4: 

£ = [. + .:^ 5 -. (13) 

* JO 


For supersonic speeds, the solution is computed using the 
Rayleigh pitot equation that is derived from adiabatic nor- 
40 mal shock wave relationships. For "pi. 4: 


and 

a {y *)=cos a cos Jq (>)t) +sin a sin X (yM co$ 


Because equation (8) is a quadratic equation, there are two 
solutions. In general, the correct estimation can be computed 
by picking the value of p that is nearest to zero. However, 
there are times when this will not work, and the full 
mathematical procedure for picking the correct root is 
detailed below in the “Mathematical Description of the 
Invention" section. 

In a preferred embodiment, the output angle of sideslip is 
calculated as the mean of the values computed using the four 
individual triples. Clearly, if one of the ports along the 
vertical meridian is deemed unusable and is weighted out of 
the algorithm, then only one valid triple remains for com- 
puting the output angle of sideslip, 
d. Airdata Module 

Once the values of the output angle of attack and the 
output angle of sideslip (a t . and p^> have been computed, 
then the incidence angles at all of the pressure ports can be 
evaluated, and only the calibration parameter 6, the 
freestreaui static pressure p M , and the impact pressure q,. 65 
remain as unknowns. However, c is implicitly a function of 
p^and q, and the resulting system of equations is non-linear. 


q t 166.92 Ml ^ (14) 

A* = [7Mi - l] 2 5 

Equation (14) is solved using a Taylor's series expansion 
and a reversion of series to obtain Mach number. 

For high Mach numbers, a very significant temperature 
rise occurs across the bow shock wave and the fluid no 
50 longer behaves as a calorically perfect gas. The large tem- 
perature rise across the shock wave causes vibrational 
modes of the air molecules to excite, which draws energy 
away from the flow. Thus, the stagnation temperatures 
actually encountered are considerably lower than those 
55 which would be computed if perfect gas calculations with 
y= 1 .4 were used. 

As a result of this non-adiabatic flow behavior, some error 
is expected to occur in the Mach number calculation using 
equation (14). Fortunately, because pressure is a mechanical 
60 quantity and depends primarily upon the mechanical aspects 
of the flow, influences caused by high temperature gas 
properties are secondary. Numerical analyses of equation 
(14) have shown that the error introduced is less than 0.2% 
for Mach numbers below 4.0. 

The airdata computed as described above by the various 
modules of the airdata estimation and evaluation system 140 
are then sent to the airdata processor 145. The airdata 
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processor 145 uses the computed airdata to determine 
whether adjustments to physical parameters of the flight 
vehicle 110 are necessary. Examples of physical aspects of 
the flight vehicle are its velocity, attitude, altitude, and 
deflection of the control and flap surfaces. In addition, the 
airdata processor 145 may send the computed airdata to the 
instrumentation or navigational systems of the flight vehicle 
110 before being sent to the control system 150 for adjust- 
ment of these flight parameters. 

111. Functional Overview of the Invention 

The airdata estimation and evaluation method of the 
present invention may be implemented in the airdata esti- 
mation and evaluation system 140 as a dedicated logic 
circuit or a field programmable gate array. Preferably, as 
shown in FIG. 1, the airdata system 140 may contain a 
microprocessor 142 and a memory 144 for storing and 
executing software coded with this method. 

The memory 144 may be a hard disk drive, CD-ROM, 
random access memory (RAM) or any other computer 
memory device or medium capable of storing the software 
or instructions containing the method of the present inven- 
tion. The microprocessor executes the method stored in the 
memory 144. 

The airdata method of the present invention is capable of 
interfacing with the FADS system 100 by receiving pressure 
measurements from the surface pressure sensor matrix 120, 
computing necessary airdata, and sending the computed 
airdata to the airdata processor 145. The airdata processor 
145 then uses the computed airdata to effect the necessary 
physical changes to the flight vehicle 110. 

Similar to the airdata estimation and evaluation system 
140, the airdata estimation and evaluation method includes 
a triples formulation function, an angle of attack function, an 
angle of sideslip function and an airdata function. 

The mathematical details of the equations given thus far 
and below are developed more fully in the section below 
entitled. “Mathematical Description of the Invention.” 

1. Triples Formulation Function 

FIG. 6 is a detailed functional flow diagram of the method 
that is preferably contained in the triples formulation module 
510. TTie method transforms a flow model equation into a 
triples formulation equation using pressure quantities from 
the surface pressure sensor matrix 120. 

The function starts with a given flow model equation. This 
equation is given by: 

sin 2 (0|)}+P^ (15) 

Next, the function determines the triples, which are defined 

as: 

r, ti *prp*' r < 1 6 > 

After the triples have been computed, the triples formu- 
lation equation is calculated. Substituting and collecting 
terms yields the triples formulation equation: 

r .kitij+b, tan tan |3] 2 +r*,(tt,+£« tan pl 2 =0 (17) 

Equation ( 17) is significant because the equation is in terms 
of the measured pressures from the surface pressure sensor 
matrix 120 only. 

The triple formulation function then sends the equation 
( 17) to the angle of attack function and the angle of sideslip 
function. 

2. Angle of Attack Function 

FIG. 7 is a detailed functional flow diagram of the angle 
of attack function of the airdata estimation and evaluation 
system of FIGS. 1 and 5. The function begins by receiving 


the surface pressure distribution from the surface pressure 
sensor matrix 120. The function then selects pressure 
quantities, p ( , p y , p k , from pressure ports lying on the vertical 
meridian 220, as shown in FIG. 2. 

5 Meanwhile, the function receives the triples formulation 
equation, equation (17). from the triples formulation func- 
tion and transforms this equation into the angle of attack 
equations. These angle of attack equations are given as: 

10 I ,M| (l8) 

<*= ■= 


and. 



The function then determines whether the magnitude of 
20 angle of attack, a, is less than or equal 45 degrees. If yes, and 
|a|^45°, equation (18) picks the correct root. However, if 
no, and |a|>45 0 , then the correct root is given by equation 
(19). These computations of the angle of attack are per- 
formed for each of the pressure ports selected by the 
25 function as previously described. 

In a preferred embodiment, the function next computes 
the mean angle of attack using the individually computed 
angles of attack. This mean angle of attack is then outputted 
to the angle of sideslip function. 

30 3. Angle of Sideslip Function 

FIG. 8 is a detailed functional flow diagram of the angle 
of sideslip function of the airdata estimation and evaluation 
system of FIGS. 1 and 5. The function starts by receiving the 
pressure distribution from the surface pressure sensor matrix 
35 120. Next, the pressure quantities, p y , p ; , p k , are selected 
from the pressure ports such that one of the pressure 
quantities is off the vertical meridian 220. 

Meanwhile, the function receives the triples formulation 
equation, equation (17), from the triples formulation func- 
40 tion. This equation is transformed by the function into the 
angle of sideslip equations, given by the quadratic equation: 

A' ian 2 ^f25' tan ^-C=0 (20) 

45 Next, the computed angle of attack is received from the 
angle of attack function. Because equation (22) is a qua- 
dratic equation, two values of the angle of sideslip, p! and 
P-, are computed from equation (22). If the magnitude of p, 
is less than the magnitude of p 2 , then the solution is equal 
50 to P,. Otherwise, the solution is equal to p 2 . The value of the 
° angle of sideslip is computed for all pressure values selected 
from the pressure distribution as described above. 

In a preferred embodiment, the computed angle of side- 
slip is calculated as the mean of all the computed values of 
55 the angle of sideslip for each pressure port. This computed 
mean angle of sideslip is then outputted to the airdata 
function. 

4. Airdata Function 

FIG. 9 is a detailed functional flow diagram of the airdata 
function of the airdata estimation and evaluation system of 
FIGS. 1 and 5. The function starts with the given flow model 
equation, equation ( 1 ): 

lcos 2 (8/)+€ (I) 

65 Once the values of angle of attack and angle of sideslip 
have been calculated, then the incidence angles at all of the 
ports can be evaluated, and only €, p^, and q t remain as 
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unknowns in ihe pressure equations. Assuming that n pres- 
sure measurements are to be used, the flow mode! equation, 
equation ( I ), can be written in matrix form with static and 
impact pressure broken-out as variables: 


COS'(fl| ) + 

a's l 




COS : (Mi) + 
esin : (tf„) 


Hr 

P,. 


( 21 ) 


ID 


Next.~the function solves for the static pressure, q r , and 
the impact pressure, p^ using weighted linear regression. 55 
The result is: 


16 

tern is that airdata states are related to the surface pressures 
by a flow model equation that captures salient features of the 
flow and is valid over a large Mach number range. The 
resulting flow model equation takes the form: 


(cos J (0 H-« 


(A-l) 


where 


cos( 0, )=cosi a, k:os( (3,. )cosi X, H-sin{ )sin< 0, )sin< X, H- 

sinla^teosl^Kos^JsiniX,) (A-2) 

Equations (A- 1) and (A-2) are the same as equations 1 and 
2 as described above. 

Defining the terms: 

a y =cos a, cos X,+&in a, sin X, cos 0,. 6=sin X, sin 0, ( A-3 ) 


Z « ~Z ,? ' (a) 

1 * 1 -t 

.-1 

Yj qM) 2 

1*1 t*l 



I q, Z <?/(U) 2 - £ qM) 

,i«i Jiw =- 1 


where £},-=cos 2 (9 ; )+e sin 2 ^). 

Meanwhile, the pressure distribution is received from the 
surface pressure sensor matrix 120. In addition, the com- 
puted angle of attack and the computed angle of sideslip are 
received from the angle of attack function and the angle of 30 
sideslip function, respectively. 

Next, the function writes equation (22) as an iterative 
estimator, given as: 


(22) equation A-2 reduces to 

cost 0, *=<3, cos((3,)+b ( sim(3,) (A-4) 

20 

a. The Triples Formulation 

By using equation (A-l) to take strategic combinations of 
three pressures( ‘‘triples"), the result is that the impact 
pressure, static pressure, and calibration parameter are 
- 5 decoupled from the angle of attack aand angle of sideslip (3. 

pi - Pj <7 c |cos 2 6, +■ esurtf, ) + - (<7<-|cos : 6X + eswr&j] + /*»] (A-5) 

Pj ~ Pk qd cos : 0, + fsin 2 # r } + - [qd cos ~$ k + esurtf* } +• Z 5 *,] 

{(1 -£)COS : 0; +£} - j(l -£)COS Z 9j + fi} 

|( 1 - £)cos- 0 y + £) - {( 1 - £)cos : 6 k + s) 

cos : 0, -COS'^j, 
eos 2 f) - cos 2 0 k 


kc 

= [[MjflQMu-X' MjfiQ) 

Pi 


Ij* II 

. Pn 


From equation (23) the function computes 6 iteratively 40 
using the values for Mach number resulting from the pre- 
vious iteration for p«, and q r . 

If the Mach number is estimated to be less than one, for 
air ("pi. 4) the Mach number is computed using the subsonic 
equation from above: 45 

< 13 > 


For supersonic speeds, the solution is computed using the 5 0 
Rayleigh pitot equation that is derived from adiabatic nor- 
mal shock wave relationships. For pi. 4: 

q c 166.92 Ml { 04) 

P*, ~ [7 M 2 , - 1 p 55 

The function solves equation ( 14) for Mach number using a 
Taylor’s series expansion and a reversion of series. 

Using these computed airdata parameters, other airdata 60 
may be calculated as needed by the FADS system 100. This 
computed airdata is then outputted to the airdata processor 
145 that determines it response as described above. 

IV. Mathematical Description of the Invention 

Following are the mathematical details of the apparatus 65 
and method of the present invention. As discussed above, the 
basic concept of the airdata estimation and evaluation sys- 


Expanding equation (A-5) and defining 

ru*prp* T pprPi' r '*FP*rPj {A ' 6) 

equation (A-5) reduces to 

cos 2 0 y +IV cos 2 0 t +r A , cos 2 0 ,O (A- 7 ) 

Substituting in from equation (A-4), dividing by cos 2 p 
(P*±90°) and collecting terms, equation (A-7) becomes 

T ik [aj+bj tan f3] 1 +Y Jl {a k +b k tan tan (3 1^=0 (A-8) 

It is important to note that in equation (A-8) the angle of 
attack a and the angle of sideslip p have been written as a 
function of the measured pressure only. Significantly, the 
impact pressure q r , the calibration parameter €, and static 
pressure P^, have been decoupled, 
b. Angle of Attack Solution 

Angle of attack is decoupled from angle of sideslip using 
a “meridian solution", where, along the vertical meridian 
<(>=0 o and <|h-± 180 o . Equation (A-8) thereby reduces to: 

r iA lcos a cos Xy+sin a sin X, cos 0 y l 2 +r y ,(cos a cos X*+sin a sin 
cos $ t | 2 +r\Jcos a cos X,+sin a sin X, cos 0,] 2 =O ( A-9) 

Factoring cos a (a*±90°) out of equation (A-9), 

r, A |cos X r +-tan a sin X, cos typ+rjcos X A -nan a sin X* cos <fcj*+ 
r\ y [cos X ( +un a sin X^ cos 0J-=O (A- 10) 

Expanding the squares, the result is a quadratic expression 
in tan a, 

[T, k sin 'X,+r,, sin*X t +r Ay sin 2 Xjlan 2 a+ 
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:ir„ cos X, sin X y cos 0,+r /V cos k k sin X A cos 0 fc +P v cos X, sin X t 
cos 0,|tan ot+ 


rr <4 COS' V r /( COS 2 X k +r k/ cos 2 A/H) 


< A- 1 1 > 


Equation {A- 14) can he reduced further by noting that 
cos 2 X=l-sin 2 X and 

Y fk cos 2 Xj+Y Jt cos 2 X A +r*, cos 2 X,— X y HT, ( i 1-sin 2 X*H- 
r i; ( I -sin 2 A.,) 

=r ik +r„+r h -< r, t sin 2 X y +r y/ sin 2 X k +r ki sin 2 X,) 

=< p,-pi H4 PrP, W PtrPj ^ L* sin 2 X y +r y , sin 2 X k +r kj sin 2 X,) 

Hr a . siirXy+F,, sin 2 X^+r^ sin 2 X,) (A- 12) 

Substituting Equation (A- 12) into Equation (A- 11), 

[T^ sin 2 Xj+Tjj sin 2 X t +r Ay sin 2 Xjftan 2 a-l[+ 

2[T lk cos Xj sin X J cos 0,+P,, cos X k sin X k cos 0*+!^ cos X, sin X, 
cos 0‘}tan ct=0 (A- 13) 

Defining, 

A=(Y tk sin 2 Xf+Yjj sin 2 Xt+r^ sin 2 X<) 

5=(T 4i cos 0, sin X, cos X y +r y , cos 0 t sin X k cot ^+Y kj cos{0 4 sin X, 

cos X 4 ) (A- 14) 


r-ir^Z+r^r+r^ 2 } <aiu 

Equation A-20 reduces to a quadratic equation in tan p 

A' tan 2 $+2B' tan (3+C=fl (A-22) 

The two solutions to equation A-22 are 


tan0|, 


& 

1 1 b ' f 

c 

<a-2: 

' *r\ 


A 7 



10 


d. Root Selection Criteria for Angle of Sideslip 
The relationship of tan[p, ] to tan[p 2 ] is developed by 
taking the product of the roots: 


15 


tan£ttan& = 


(A-24) 




C 

' A 7 
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Rearranging equation (A-24) gives the relationship of the 
two roots 


Cl C r* , 

25 “7“l2" fc l 


(A-25) 


Equation (A- 14) is written in the simple form 
A [tan 2 a-l]+2i? tan ct=0 


(A- 13) 


The “correct” solution is extracted from equation (A- 15) by 
noting that: 


tan - a - 1 = 


sin 2 a - cos 2 a 1 - 2cos~a cos2a 


(A-16) 


cos-a cos-ar 


For several specific port arrangements relevant to the 
surface pressure sensor matrix design as described 
30 previously, the relationship between the two roots can be 
written in closed form. Three of these special arrangements 
will be now be discussed. 

e. Root Selection Criteria for the Triples on Lateral Meridian 
Along the lateral angle-of-sideslip meridian, the outer 
35 port clock angles are ±90°, and a J>y .=cos a cos X iJk , b ijk =±sin 
X iJk . , thus equation (A-24) reduces to 


Substituting equation (A- 16) into equation (A- 15), and 
multiplying through by cos 2 a. 


-A cos 2ct+25 cos a sin a=-A cos 2a +5 sin 2a^=0 

and the solution for a is: 


1 -iMl 
a= 5 un [-j 


(A- 17) 


(A-18) 


40 


[Cl r,ta 2 + f>fl£ + r k ,aj 
[a 7 ]" r lk bj + r Jt bi + v kJ trr 

, r^COS 2 A + r ; ;COS 2 A* + r* y COS 2 A # 

= cos" a — — — — —77 r — “v, 

T.jtsin-A, + r y ;sin*A* + r* ; surA, 

= -cos 2 a 


(A-26) 


45 


For |a|^45°, equation (A-18) picks the correct root; when 
|ot|>45 0 then the correct root is given by the complement of 
equation (A-18): 


Thus, along the lateral meridian, the relationship between 
the two solutions of equation (A-22) is: 


50 


tan/?t = cos 2 tttanj$: - - J 


(A- 27) 




(A- 19) 


c. Angle of Sideslip Solution 

Given the solution for angle of attack, a procedure similar 
to that of the previous section is followed for angle of 
sideslip. Expanding the squares in equation (A-7) and solv- 
ing for like terms in tanp, 

[r>/+r,A 2 +r*A 2 |tan 2 (3+21 Y lk u / b ) +Y Ji a k b k L'+Y kJ a i b i I tan 

^r ik xi'+Y ){ a k ^Y k/ ar=0 ( A-20) 


Even for angles of attack as high as 50 degrees the two 
angle of sideslip solutions are nearly 70 degrees out of 
phase. Consequently, for triples along the lateral meridian, 
55 selecting the root with |^| closest to zero will always gives 
the correct solution for (5. 
f. Off-Meridian Triples, Symmetric Arrangement 
Assume ports i and j lie the lateral meridian and port k lies 
on the vertical meridian, then 

60 


Defining 

a= ( i 1 V>/+r ; A 2 +r I 

S' 3 ( Y^fa+Y^b^Y^b, ) 


65 


fa t = ^ , => a ; = cosacosA,, b; = sint, 
ss - , =* a } = cosmrosA, , b f - - sinA } 

-fljtscos a cos X*+sin a sin X^ b k =0 


(A-28) 


(A-28) 
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and for a symmetric arrangement, X—X,. and C'/A’ reduces Factoring out cos" a 


to 


c , | LosA“(r* > r i; M- r /t (cosA* f tajkr.iaV | (A 29) 

— = cos') f c ; — 

v \ sm-A/ir,* + r*,j 


Noting that r^+r^p-r,, - Equation (A-29) can be reduced 
further 


cos-(ir-U) 

c* r.i + r ; — + r^cos-i 

c 1 1 COV-U 

— = cc»s‘a< r: 

A' I fusin' A, 


cos*(a - X k ) > 

r* > r fi — ; — ; — + r Ay < i - sin - a, ) 


(A-35) 


10 


c 

— = COS’ 1/ 1 

A' 


~r y( cosA’ + f j, (cosA* + tantrsinAt ) : 


- f H sin - A , 


(A- 30) 


cos*<r 
sin 1 


— ; — cos A* — (cosocosA* + simsinA* )’ 

iirA, L J cos- «• 


cos'ocosA' - cos’ftf - /V* ) 
siirA, 


Substituting Equation (A- 30) into equation < A-25), the rela- 
tionship between the two roots becomes 


15 


20 


Collecting terms, equation (A-35) reduces to 

cos‘(<j — A* ) 


C 

— = COS' O' 

A' 


(f iA +ri y ) + v jt - 


r*,(sin - A,) 


(A-36) 


sin’ A, 


cos ’(a - A* ) _ . , ' 

r -> + r /- - r * ( * urA> 


r i; sin : A 


tan/?i = 


cos : crcosA^ -cos : (a-A*) 


siirA, 


tan[ - - | 


(A-31) 




( (c°s^-A.)_i| j 


2 

J cos-o J 

V ' 1 


= cos a 

siirA, 

25 



For the surface pressure sensor matrix, there are two The pressure dependent term, 
possible symmetric off-meridian angle of sideslip triples 
combinations, 30 


{A„ X,, X*}={20°, 20°, 45°), and {X,, X,, X*}={20°, 20°. 20°lJA-32) 


£> 

W 


* = 


sinA; 


(A-38) 


For the first set of coordinates {X,, X yt X it }={20 0 ,20 0 ,45 0 } the can be removed from equation (A-34) by substituting in 
system is highly ill-conditioned. This ill-conditioning make 35 from equation (A-5). Defining 
selecting the correct root impossible. Furthermore, when the 

second set of coordinates {X,, X y , XJ={20 0 ,20°,20 0 } is used [cos : <a -A*> J (A-37) 

a singularity exist at a^0° and P=0°. 

Clearly, using the X*=45° port is the preferable choice for 
the off-lateral meridian. If one assumes that true angle-of- 40 
sideslip limits to be experienced are less than ±20°, then for 
17°,and a ^20°, a “root closest to zero” criterion is valid 
for selecting the “correct” root. For values of a between 1 7° 

and 20°, the X { , y(t} ={20°,20 0 ,45°} system is indeterminate, — = cos-erf — ^ - lj 

and instead the X {iyjt} ={20 o ,20 o ,20°f triple is used. With this 45 A * r * ; * 

substitution, then a “root closest to zero” criterion always 
selects the correct root. 

g. Off-Meridian Triples, Asymmetric Arrangement 

The root locations for one special case asymmetric port 
arrangement can be analyzed, that is the arrangement where 
Xy=0. Consider the case where, 

0, = ± * , ^ a t = costrcosA, , b; = ± sinA, 

$ y -=0. Ay=tHay=5eos a, b ~ 0 

0 A =O. *-« A =cos a cos X A +sin a sin X A cosUa-X*). b k = 0 (A-33) 

Following the same procedures as before 


50 


then equation (A-37) reduces to 


and finally 

tan/? | = cos“tf| ~~ - 1 |tan{ — -&] 


(A-39) 


For the surface pressure sensor matrix 120 design, there 
are four possible asymmetric off-meridian angle of sideslip 
55 triples combinations, 

{{K <M.(V MiK ^}>={{20°.±<Xn. {0 U .0°}, {45°.0°} }(A-40) 


and 


60 


UK ♦idl-darton. {o u .o°}. rso-jtnXA-m 


C r,*cos'a + r /4 cos - (a - A* ) + r^cos’u-cos-A, 
A' r\y siirA; 


For the first sets of coordinates is {X^, 4> v , <t> ; }, {X A , 
4>,}}={{20° +90°}, {0°,0°},{45 o ,0 o } }. In the range of 
|p|<25°, it is very interesting that the second solution is 
(A -34) nearly independent of the value of the first solution. For 

65 <^-=90°, p 2 2-80°; and for <(>p-90 o , p 2 s8°. 

The second set of coordinates is {X,, <(>,,}, {X y , fy}, {X*, 
4>^1 }— ( {0°,0°}, {20°,0°}}. Again in the range of 
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|p|<25°, it is very interesting to note that p 2 is identical to the 




value given by the previous set of coordinates (X*=45°). 
Clearly, the terms involving X k are canceled out of equation 


X * -t ® ux ' ) 

£ aihf.r. 

(A-35). As a result for the asymmetric configuration, a “root 


n n 


closest to zero” criterion will always pick the correct solu- 5 

[*1 


- 7 q,<ri) 5"! 


tion. 



<>i 


cos ~( 0 /) + 
cos’ (tin) + 


iO 


15 


I A-42) 


25 


30 


h. Angle of Sideslip Solution Summary 

The quadratic equation, (A-22), has two roots, with no 
guarantees that they be in general orthogonal. However, 
several conclusions can be made though for special geom- 
etry arrangements that the surface pressure sensor matrix 
120 exploits in its design. These conclusions are: 

i. For triples which lie on the lateral meridian, the solu- 
tions are nearly orthogonal, and for |p|<25°, using a 
'‘root closest to zero” selection criterion always pick the 
correct solution. 

ii. For off-meridian triples arrangements, only the sym- 
metric arrangement presents a problem. For this 


all three incidence angles are equal. In the vicinity of 
this singularity, the system is highly ill-conditioned. 
This ill-conditioning makes the system indeterminate 
and selecting the correct root is impossible without 
prior information. 

iii. For the surface pressure sensor matrix 120, using the 
lower 45° sensor location to complete the symmetric 
triple puts the singularity at 18.207° local angle of 
attack. Using the lower 200 sensor location puts the 
singularity at 0° local angle of attack. 

iv. If one assumes that true angle of sideslip limits to be 
experienced are less than ±25°, then for the ^{/y*} = {20°, 
20°, 45°} triple, for 17°, a^20°, a “root closest to 
zero” criterion is valid for selecting the “correct” root. 

For values of a between 17° and 20°, the \ {ijk} ={20 ° , 
20°, 45°} system should be considered indeterminate, 
and instead the A, {,•,*}={ 20°, 20°, 20°} triple should be 
used. If this substitution is made, then a “root closest to 
zero” criterion will always be valid for selecting the 
“correct” root. 

v. For the asymmetric off-meridian configuration, the 
alternate p solution will always hover around ±80 
degrees. Thus the “root closest to zero” criterion is 
always valid. 

i. Static Pressure, Impact Pressure, and Mach Number 
Solution 

Once the local values of a and p have been determined, 
then the incidence angles at all of the sensor locations can be 
evaluated, and only e, p M , and q c remain as unknowns in the 50 
pressure equations. Assuming that n pressure measurements 
are to be used, the flow model equation (equation ( A- 1 )) can 
be written in matrix form with static and impact pressure 
broken-out as variables: 


35 


40 


Ignoring for now the fact that E is implicitly a function of 
the static and the impact pressure, equation (A^t2) can be 65 
solved for the static and impact pressure using weighted 
linear regression. The result is: 




(A-43) 


I m I - I *(«.) 


where Q=cos 2 (0,)+€ sin 2 (0 4 ). In equation (A-43), the q/s 
are weights which have a nominal value of 1.0. Setting the 
value of q, to zero weights the i'" 1 pressure reading out of the 
algorithm. Equation (A-43) looks like a closed form 
solution, but the € terms on the right hand side of equation 
( A^O) are implicitly a function of p^ and q c . Thus equation 
( A-43) is actually non-linear and the solutions for p M and q c 
must be extracted iteratively. Defining 


(A-44) 


where Q / </) =cos 2 (0 / >+€ </) sin 2 (0 4 ), equation (A-43) is written 
as an iterative estimator 


20 

(nV 1 ) i 


<l\ 

- 0 

M<j) = 

(H^) I 


0 - 

•' <Jn 




P l 


(A-45) 


In equation (A^f5) the subscript (j) refers the result of the 
j rh iteration. The matrix M w is defined in equation (A-44) 
with e being evaluated using the values for Mach number 
resulting from the previous iteration for p„ and q c . 
j. Algorithm Stability Analysis 

Since equation (A-45) is to be implemented as part of a 
real-time airdata estimation algorithm, it is essential that 
potential instability regions be identified. In general, it is 
impossible to analyze the stability of non-linear equations in 
two unknowns. However, a linearized stability analysis will 
determine the behavior of the system with respect to small 
disturbances. Defining the terms: 


45 


P i 


(A-46) 


and recalling that € is a function of Mach number, and 
implicitly a function of static and impact pressure, then 


55 


M[x\ - 


n, i 




(A-47) 


60 


and equation (A-42) can be re-written as; 
2=M[x\k 


(A-48) 


Using the definitions of equation (A-46) and (A-47), the 
iterative estimation equation can be written in the form: 

(A ^ 9 > 

In equation (A-49), is the estimate after the j+T" 1 
iteration, and M[x 0> ]is the matrix of equation (A-47) evalu- 
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ated using the result from the jth iteration. Subtracting 
equation (A-49) from (A-48). and expanding M[x] in a 
Taylor* s series about and neglecting terms higher than 
first order in the perturbations, equation (A-49) reduces to: 

S 

( , 7Af y ' ') »}( jt— 1 *>=0 ( A-50) 

Defining the error vectors, 

x^' '*=i v-x ^‘ ') ( A-5 h 

10 

Equation (A-51) becomes the linearized error equation for 
the iteration 

’=0 < A-52) 

15 

and the eigenvalues (roots) of the characteristic equation 

Det{ ([Af[.v v> ] r Ml-^’]r 1 ( A-53 ) 

will determine the linear-stability of system. The elements of ^ 
the Jacobian, VM°\x, are evaluated as 

f . - d c (J) . d'(f) d\i M 1 (A- 54) 

I 9cSm-9, — — q c sm-fl, — — | 


' <3 AC dq c 


. . d'(J) dM„ . , d f </> dM x 

* s “ rS ' JmZ JK. 


In equation (A-54) the parameter 3Q 

A. 

DM m 

is the sensitivity of the calibration parameter to Mach 35 
number and is determined by numerically differentiating the 
calibration parameter (eM) with respect to Mach number. 
The parameters 

DM„ 40 

dq c 


and 


3M m 

dP M 


45 


flare evaluated by differentiating the isentropic flow equa- 
tion for subsonic flow, or the Rayleigh Pitot equation for 50 
supersonic flow. For y^l.4, the resulting expressions are 


DM l DM 

J^ c "~W'dK, 

M DM 



(A-55) 

55 


where the Mach number derivative is 


(— ) 

\DM / iuitfOlUC 

= 1.4 AC[I +0.2 Mi \ : 

tor M» 

s j (A-56) 60 

Df \ 

. 12 M; 

, -I] , 


— — J 

= 1 168.45 Mi 

: -7 for 

M m > 1 

DM f tuf Mrwruc 

17 Mi 

- *1 5 



65 

Substituting equation (A-55) into equation (A-54), the Jaco- 
bian matrix reduces to 



Equation (A-57) is dependent only on the freestream Mach 
number and the local incidence angles. Therefore, for a 
given surface pressure sensor matrix geometry, the stability 
of the system (eigenvalues equation (A-53)) is dependent 
only on the freestream Mach number and the local flow 
incidence angles. This conclusion is very important, as it 
allows the stability characteristics of the system to be 
analyzed independent of trajectory. 

In order to achieve the landing and ascent airdata require- 
ments for the X-33 flight vehicle, the airdata estimation and 
evaluation system of the present invention is capable of 
operating between Mach 0.20 and Mach 4.0 and meets the 
following accuracy requirements: 

i) Mach Number: ±5.0% accuracy for 2.5 ^M ^4.0, 
±2.50% accuracy for 0.6 = M<2.5, ±0.015 absolute 
error between M 0.20 and M 0.60; 

ii) Angle of attack: ±1.5° absolute accuracy for first 3 
flights, ±0.50°, thereafter; 

iii) Angle of sideslip: ±0,5° absolute accuracy for all 
flights; 

iv) Geopotennal altitude: ±200 feet absolute accuracy for 
0.2^M^4.0; and, 

v) Dynamic Pressure: ±15 psf. for 0.2^ M ^4.0. 

V. Conclusion The airdata estimation and evaluation system 
of the present invention computes the airdata and flow 
parameters from nonintrusive surface pressure measure- 
ments. This system and method transforms a flow model 
equation into a triples formulation equation whereby airdata 
may be extracted without the use of nonlinear regression 
methods. This results in a more accurate and more stable 
system and method of estimating and evaluating airdata in a 
FADS system. Moreover, the code required to implement 
and maintain the method of the present invention is less 
complex than prior art algorithms. 

The foregoing description of the preferred embodiment of 
the invention has been presented for the purposes of illus- 
tration and description. It is not intended to be exhaustive or 
to limit the invention to the precise form disclosed. Many 
modifications and variations are possible in light of the 
above teaching. It is intended that the scope of the invention 
be limited not by this detailed description, but rather by the 
claims appended hereto. 

What is claimed is: 

1. A flush airdata sensing (FADS) system, for use on a 
vehicle in a flowfield, comprising: 

surface mounted pressure sensors located on the vehicle 
for providing a surface pressure distribution; 

a airdata estimation system, connected to the surface 
mounted pressure sensors, for estimating and evaluat- 
ing airdata from the surface pressure distribution and 
the flow model equation, the airdata estimation system 
comprising: 

a triples formulation module that transforms the flow 
model equation into a triples formulation equation 
by using pressure quantities corresponding to the 
surface pressure distribution; 
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an angle of attack module that computes an angle of 
attack quantity from the triples formulation equation 
by choosing the pressure quantities such that the 
angle of attack quantity is decoupled from an angle 
of sideslip quantity; 

an angle of sideslip module that computes the angle of 
sideslip quantity from the triples formulation equa- 
tion by using the angle of attack quantity; 
a airdata module that estimates and evaluates airdata by 
using the angle of attack quantity and the angle of 
sideslip quantity; and, 

a control system computer, connected to the airdata 
estimation system, for adjusting an orientation of the 
vehicle in the flowfieid from estimated and evaluated 
airdata. 

2. The invention as set forth in claim 1, wherein the flow 
model equation is a splice of a subsonic potential flow 
solution for a blunt body and a hypersonic modified New- 
tonian flow. 

3. The invention as set forth in claim 1, wherein the triples 
formulation module further comprises: 

a triples difference module that transforms the flow model 
equation into the triples formulation equation by taking 
the difference of pressure quantities from the surface 
pressure distribution. 

4. The invention as set forth in claim 3, wherein the 
pressure quantities are three pressure quantities from the 
surface pressure distribution. 

5. The invention as set forth in claim 3, wherein: 

the flow model equation is a function of a first pressure 
quantity of the surface pressure distribution, a static 
pressure, an impact pressure, a position error calibra- 
tion parameter, the angle of attack quantity and the 
angle of sideslip quantity; and, 

the triples formulation equation is a function of the first 
pressure quantity, a second pressure quantity, and a 
third pressure quantity, all corresponding to the surface 
pressure distribution, and the angle of attack quantity 
and the angle of sideslip quantity. 

6. The invention as set forth in claim 1, wherein the angle 
of attack module further comprises: 

a pressure selection module that chooses the pressure 
quantities such that they ail lie in a vertical plane of the 
surface pressure distribution. 


7. The invention as set forth in claim 6, wherein the 
number of pressure quantities is at least three. 

8. The invention as set forth in claim 1, wherein the angle 
of attack module further comprises: 

s a mean AOA module that computes a mean angle of attack 
quantity from each computed angle of attack quantity 
of the surface pressure distribution. 

9. The invention as set forth in claim 1, wherein the angle 
of sideslip module further comprises: 

10 an AOS equation module that transforms the triples 
formulation equation into an angle of sideslip equation. 

10. The invention as set forth in claim 9, wherein the 
angle of sideslip equation is a quadratic equation. 

11. The invention as set forth in claim 10, wherein the 

15 quadratic equation comprises a trigonometric function of the 

angle of sideslip quantity. 

12. The invention as set forth in claim 11, wherein the 
trigonometric function is a tangent. 

13. The invention as set forth in claim 10, wherein the 

20 angle of sideslip quantity is one of two roots of the quadratic 

equation. 

14. The invention as set forth in claim 13, wherein the 
angle of sideslip quantity is a root of the quadratic equation 
having an absolute value of the root closest to zero. 

25 15. The invention as set forth in claim 1, wherein the 

angle of sideslip module further comprises: 

a mean AOS module that computes a mean angle of 
sideslip quantity from each computed angle of sideslip 
quantity corresponding to a pressure of the surface 
pressure distribution. 

16. The invention as set forth in claim 1, wherein the 
airdata module further comprises: 

a matrix module that writes the flow model equation as a 

35 matrix to form a flow matrix equation and solves the 
flow matrix equation for a static pressure and an impact 
pressure. 

17. The invention as set forth in claim 16, wherein the 
matrix module further comprises: 

40 an iterative estimator module that constructs an iterative 
estimator from the flow matrix equation and solves 
iteratively for the static pressure and the impact pres- 
sure. 


* * * * * 


